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1 Liebmann technical documentation series 


1. Wyznaczanie rozktadu pola elektrostatycznego w prézni metodg relak- 
sacyjna Liebmanna. (Polish version / wersja polska) 


2. Determination of electrostatic field distribution by using Liebmann relax- 
ation method. (English version / wersja angielska) 


3. Graphics. Mapping voltages to colours (colormaps). 


4. Laplace equation 2D (XY). (Cartesian coordinates). Relaxation scheme 
explained. (5 - point star) 


5. Laplace equation 2D (ZR) (Cylindrical coordinates). Relaxation scheme 
explained. (5 - point star) 


6. Liebmann source sode. (ANSI C programming language) 


2 Versions of this document 


1. version 1 - 2023.11.03 
. version 2 - 2024.01.26 
. version 3 - 2024.02.02 
. version 4 - 2024.02.05 
. version 5 - 2024.05.18 
. version 6 - 2024.05.23 


N ODO Oo FF WO DPD 


. version 7 - 2024.05.24 


3 Solving Laplace equation using relaxation method 


| tried to solve Laplace equation using mainly information from Pierre Grivet’s 
book (Electron Optics) - [1]. 
There are few editions of this book (1965, 1972). Second edition (1972) con- 
tains explanation of relaxation method (page 38). 

More generalized approaches has been drafted by James R. Nagel - [2]. 


https://my.ece.utah.edu/~ ece6340/LECTURES/Feb1/ (visited 2023-03-01). 


There are also publications edited by Albert Septier: Focusing of Charged 
Particles [3] and Applied Charged Particle Optics (part A). [4]. 
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| have also found some ideas in publication of D W O Heddle: Electrostatic 
Lens Systems [5] (especially using PC computers to solve electrostatic prob- 
lems). 

| have also found (brief) description of by - hand solving of Laplace equa- 
tion by Bohdan Paszkowski - [6] (Polish edition). English translation of this book 
also exists - [7]. 


| would like to thank many people, who helped me with this challenge. Espe- 
cially prof. dr hab. Mieczystaw Jatochowski (Supervisor of my master’s thesis), 
who enabled me to use SIMION and MATLAB software while writing master’s 
thesis about electron optical systems at University of Maria Curie - Sktodowska 
in Lublin in 2008. | would also thank to prof. Marcin Turek for fruitful discus- 
sion about numerical methods. What is more, my colleague Bartosz in 2012 
had explained me general problems with software efficiency. So he had also 
contributed significantly to the idea of Liebmann software (especially using C 
language). 


4 Explanation of symbols in calculations 


¢ P, - i-th mesh node 

+ V; - value of electrostatic potential at node P;. Unit - [V] 

* h- mesh step (for example h,, - mesh step in x direction). Unit - [mm] 
* gi+/— > gradient in direction i (for example giz— = ae . Unit - [~~] 


* irow - index of row in mesh. Values of 7,-¢,, = 1, 2, .., size_row 


* ico - index of column in mesh. Values of 7..; = 1, 2, .., size_col 


11 9 Mesh XY -typeA 
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Figure 1: Mesh XY type A 
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6 Mesh XY - type B 


hy F hy 
gradient V outside a mesh does not exist 
A 
y trow 
P; Px Po 
Ymax + 3 6 @ @ 
Py Ps Ps 
+2 ® © ) 
hy 
P, Py P3 
Ymin + 1 a @ @ 
hy 
1 2 3 teol 
Lmin Tmax x 


Figure 2: Mesh XY type B 


iz 7 Mesh XY - type C 
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Figure 3: Mesh XY type C 


io & Mesh XY - type D 


11 hy =hy=h 
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Figure 4: Mesh XY type D 


163 


164 


165 


166 


167 


168 


169 


170 


171 


172 


173 


174 


175 


9 Example of A-type mesh in ANSI C 


Example of A- type mesh in ANSI C program. The mesh is represented by 2 
dimensional array of double precision numbers. Rows and columns in mesh 
are numbered from 1 (this was my choice) instead of default 0 (as usual in C 
language). This choice nas pros and cons. Is is easier to calculate mesh size 
(size_row * size_col). Access to each node can be also more intuitive, but logic 
in each library function must contain this shift between node ordering styles. 


— +| 0 1 2 3 4 5 
Dy+ 
trow 
4 3 18 19 20 21 22 23 3 
3 2 12 13 14 15 16 17 2 
2 1 6 i 8 9 10 11 1 
1 0 F 0 1 2 3 4 8) 0 
IL -¥e 
Gu Gat 
—| 0 1 2 3 4 5 
Gy- 


1 2 3 4 5 6 tol 


Figure 5: ANSI C - mesh XY type A 


* gx— = double* ptr_gX_minus 
* 9r+ = double* ptr_gX_plus 
* gy— = double* ptr_gY_minus 


* gy+ = double* ptr_gY_plus 


V = double* ptr_V 


unsigned int size_row == 
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unsigned int size_col == 


unsigned int i_row == 1, 2, .., 4 


unsigned int i_col == 1,2, .., 6 


double h_x == 1.0 [mn] 


double h_y == 2.0 [mm] 


The following picture describes analogous version of ptr_V mesh, which 
can be dynamically allocated on heap by pointer metod. The mesh is rep- 
resented by single block of memory. The numbers or rows and columns are 
also known, so each node can be also accessed by appropriate index (memory 
address). 


ptr_V——+ | 0 1 2 Pe 23 


Each mesh point has its unique index (let’s say icp - (index of central 
point)), which can be determined, if we know indices of row and column (i_row, 
i_col). 


icp == ( i_row - 1) * size_col + i_col - 1 (9.1) 


For example for each point of a mesh indices of row and column have val- 
ues: 


l_row == 1, 2, .. , size_row 


i_col == 1, 2, .. , size_col 
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Example of B- tyoe mesh in ANSI C program. The mesh is analogous to A - 


Example of B-type mesh in ANSI C 


type mesh. There are no electric field gradients on mesh borders. 


trow 
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Figure 6: ANSI C - mesh XY type B 


V = double* ptr_V 

unsigned int size_row == 4 
unsigned int size_col == 6 
unsigned int i_row == 1, 2, .., 4 
unsigned int i_col == 1,2, .., 6 
double h_x == 1.0 [mn] 


double h_y == 2.0 [mm] 
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22 11 Example of C-type mesh in ANSI C 


203 Example of C- type mesh in ANSI C program. The mesh is analogous to A - 
204 type mesh. Just mesh mesh step hz = hy = h. 


—)| 0 1 2 3 4 5 
Gy+ 
brow 
4 3 18 19 20 21 22 23 3 
3 2 12 13 14 15 16 17 2 
2 1 6 7 8 9 10 11 1 
1 0 / 0 it 2 3 4 5 0 
i: 
In Ja+ 
— | 0 1 2 3 4 5 
Iy— 
1 2 3 4 5 6 tol 
Figure 7: ANSI C - mesh XY type C 
205 * gr— = double* ptr_gX_minus 
206 * 9r+ = double* ptr_gX_plus 
207 * Qy- =double* ptr_gY_minus 
208 * Gyt+ = double* ptr_gY_plus 
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V = double* ptr_V 
unsigned int size_row 
unsigned int size_col 


unsigned int i_row == 
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* unsigned int i_col == 1,2, 


* double h == 1.0 [mn] 
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2s 12 Example of D-type mesh in ANSI C 


21s Example of D- type mesh in ANSI C program. The mesh is analogous to B - 
217 type mesh. Just hy = hy = h. 


trow 

4 18 19 20 21 22 23 
3 12 13 14 15 16 17 
2 6 7 8 9 10 11 
1 / 0 i 2 3 4 5 


1 2 3 4 5 6 tol 


Figure 8: ANSI C - mesh XY type D 


218 ¢ V = double* ptr_V 

219 * unsigned int size_row == 4 

220 * unsigned int size_col == 6 

221 * unsigned int i_row == 1, 2, .., 4 
220 * unsigned int i_col == 1,2, .., 6 


double h == 1.0 [mn] 
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13 Relaxation formula for node P1 


13.1 Node description 


Left, botton corner of mesh XY. 


13.2 Calculation of relaxation formula 


Laplace equation at node P; 


Vv? (View) p, =0 


OV» Vi» 
( a2) +( “| = 6 
xv Oy 
Py Pi 


d 
| 


Let us substitute approximations to Laplace equation. 


Va-V giz  Va-Vi Qty 


| — 0) 
lige hy hy? hy 


Let us find Vy 
VY, =? 


Y-V Va-Vi Qa , Jly- 


fe. digs aig = ny 


Let us multiply both sides 


272 
| gh, 
We obtain 


Vohi, — Vihy + Vani — Vibe = gie—hehiy + giy—heh! 


Vi (he + hg) = Vahy + Vahe — gie—hahy, — Gy—hzhy 
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(13.1) 


(13.2) 


(13.3) 


(13.4) 


(13.5) 


(13.6) 


(13.7) 


(13.8) 


(13.9) 


(13.10) 
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13.3 Final forms of relaxation formula 


13.3.1. xyLV_RELAX5 P1 A 


hy # hy 
Jla—s Gly— Fx 0 
ee Vahe + Vah2 _ Jia—hehi, _ Gig2nehy 
' hz + he 
13.3.2 xyLV_RELAX5 P1B 
hy & hy 
Jla—> GPly- = 0 
Voh2 + Vih2 
CS ae 
er "y 
13.3.3 xyLV_RELAX5 P1_C 
hy = hy =h 
Jla—s Gly Fx 0 
irs V2 + V4 - Qiz—h = Piy—h 
: 2 
13.3.4 xyLV_RELAX5 P1_D 
he = hy =h 
Qa->Jly- = 0 


17 


(13.11) 


(13.12) 


(13.13) 


(13.14) 
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14 Relaxation formula for node P2 


14.1 Node 


description 


Bottom edge of mesh XY. 


14.2 Calculation of relaxation formula 


Laplace equation at node P, 


2 _ 
V (Vie.y)) p, =0 (14.1) 
PV y) PV y) 
SS Ss = 0 14.2 
( Ox? Oy? ee) 
P2 P2 
Approximation of partial derivatives of V(,.,) at node P, 
PVin,y) Be Myce -_ tack 7 Vi +4 V3 = Ws ies 
Ox? ~ hie = h2 " 
P2 
V5—V: V2—Vay— 
PV aca a hy 2 Ry y 7 V5 —V2  Gay— (14.4) 
Oy? es hy he hy 
Let us substitute approximations to Laplace equation. 
—2¥, V5-—V: = 
oes 2, Vs ne G2y- _ 9 (14.5) 
hy ity hy 
Let us find V2 
Vo =? (14.6) 
Vet —2Vo.. Va—Ve_. Gays 
= (14.7) 
h2 he hy 
Let us multiply both sides 
| Aone (14.8) 
We obtain 
Vine + Vgho — 2Vahs + Veh = goyhihy (14.9) 
V2 (he +2) = (Vi + V3) ho + Veh — goy—hohy (14.10) 
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20 14.3 Final forms of relaxation formula 
21 14.3.1 xyLV_RELAX5 P2 A 
hig Filey 
Jay— FO 


(Vi oF V3) he oh Vsh2 ~ Goy—h2 hy 
AZ + he 


Y= (14.11) 


oso 14.3.2 xyLV_RELAX5 P2 B 
he # hy 


J2y- = 0 


(Vi + V3) h2 + Veh 
AZ + hd 


2= (14.12) 
23 14.3.3 xyLV_RELAX5 P2_C 
hy =hy=h 
Gay— #0 


Vi+Va+Vi—goy—h 
eee ne 92 (14.13) 


ose 14.3.4 xyLV_RELAX5 P2_D 
he = hy =h 
J2y- = 0 


_UtVBtVs 


? 3 


(14.14) 
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15 Relaxation formula for node P3 


15.1 Node description 


Right, botton corner of mesh XY. 


15.2 Calculation of relaxation formula 


Laplace equation at node P; 


Vv? (Vey) p, = 


Vow) she View) = 
Ox? Oy? 
P3 P3 


Approximation of partial derivatives of V(,,) at node P3 


Ox? - 


Viiv. ay 
PVixy) a on Vee _ Ga V2 — V3 
7 he he he 
3 


d 


aye = 


Ve—V: V3 —V3y— 
(S52) Fy : hy j _ Ve V3 gay- 
Bs hy he hy 


Let us substitute approximations to Laplace equation. 


93at+ V2 — V3 Ve — V3 J3y- 


=0 
hy he he hy 
Let us find V3 
V3 =? 
Vo=Ve., Ver Vs. 93y= 3x4 
h2 he hy he 
Let us multiply both sides 
272 
| -hih, 


We obtain 


Vahi — Vahi + Veh — Vuh = gay-hehy — g32thahs 


V3 (AZ + ho) = Vohe + Veh + gae4 hehe — gay—hohy 
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(15.1) 


(15.2) 


(15.3) 


(15.4) 


(15.5) 


(15.6) 


(15.7) 


(15.8) 


(15.9) 


(15.10) 
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15.3. Final forms of relaxation formula 


15.3.1. xyLV_RELAX5 P3 A 


hy # hy 
932+ 1 I3y— _ 0 
We Vahi + Voh2 ao 93x4+heht _ J3y—N2 hy 
he he 
15.3.2 xyLV_RELAX5 P3 B 
hy # hy 


932+ 5 J3y-— = 0 
Vohy + Vehi 


V3 = 
h2 +12 


15.3.3. xyLV_RELAX5 P3_C 
a ee 
932+ 1 I3y—- # 0 


= Vo + Ve + g3a+h — 93y—h 


ie 2 


15.3.4 xyLV_RELAX5 P3_D 
eh 
932+ 5 J3y-— = 0 


_ Ve+Ve 


V3 5 
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(15.11) 


(15.12) 


(15.13) 


(15.14) 
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16 Relaxation formula for node P4 


16.1 Node description 
Left edge of mesh XY. 


16.2 Calculation of relaxation formula 


Laplace equation at node P, 


Vv? (View) p, =0 


OV» Vi» 
cad abe ead = 
Ox Oy 
P4 P4 


Approximation of partial derivatives of V(,.,) at node P4 


wm => 


V5—V. V. —V. Ye 
(52) iat Fe V5 Vt te 
Ps 


V7-Vi Va-Vi 
OV) ied Fey oa he ; _ Vice 24 
Oy? a hy h2 


Let us substitute approximations to Laplace equation. 


Vs—Va Q4e— Vit V7 —2V4 


| =0 
ie hy ie 


Let us find V4 
V4 =? 
Vs5—-Va  VWtVe-—2V4 — Gae- 
h2 he hy 
Let us multiply both sides 


x 


242 
| -h hy 
We obtain 


Vsh2 — Vah2 + Vih2 + Veh? — 2V4h2 = gar—hah? 


Va (2h2 +h?) = (Vi + Vz) h2 + Vsh2 — gaa—hehé 
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(16.3) 


(16.4) 


(16.5) 


(16.6) 


(16.7) 


(16.8) 


(16.9) 


(16.10) 


20 16.3 Final forms of relaxation formula 


21 16.3.1 xyLV_RELAX5 P4 A 


hice 
J4a— # 0 
oe (Vi + V7) h2 + Veh? — gar—heh? (16.11) 
40> : 
2h2 + h2 
22 16.3.2 xyLV_RELAX5 P4 B 
hie Ett 
J4e— = 0 
(Vi + Vz) he + Vohg 
V2 = - (16.12) 
2h? + h2 
23 16.3.3 xyLV_RELAX5 P4_C 
hig hype h 
J4a— # 0 
VitVs+V7—ga2—-h 
aes i ep (16.13) 
2a 16.3.4 xyLV_RELAX5 P4_D 
Cah 
G4e— = 0 
Vwt+V; 
yes f (16.14) 
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17 Relaxation formula for node P5 


17.1 Node description 


Node inside a mesh XY. 


17.2 Calculation of relaxation formula 


Laplace equation at node P; 


Vv? (Vey) py = 


View) abe View) = 
Ox? Oy? 
Ps Ps 


Approximation of partial derivatives of V(,.,) at node Ps; 


(73) - Mok _ a 7 Vii =o 
Ox? - 
Ps 


Q 


he h2 


2 


Vg—V5 V5—V: 
(52) hi aly = es Vo ae Vg = 2V5 
m hy he 


Let us substitute approximations to Laplace equation. 


VatVe—2V5 | Vat Vg— 2V5 


=0 
2 2 


Let us find Vs 


V5 =? 


=0 


V4 + Ve — 2V5 se V2 + Vg — 2V5 
h2 he 
Let us multiply both sides 


272 
| hgh, 
We obtain 


Vihi + Veh2 — 2Vsho + Voh2 + Vgh2 — 2V5h2 = 0 


2Vs (hz + Ae) = (Vo + Ve) ho + (Va + Vo) ne 
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(17.1) 


(17.2) 


(17.3) 


(17.4) 


(17.5) 


(17.6) 


(17.7) 


(17.8) 


(17.9) 


(17.10) 


295 


296 


297 


298 


299 


300 


302 


303 


304 


305 


306 


17.3 Final forms of relaxation formula 
17.3.1 xyLV_RELAX5 P5_A 
hin hig 
No gradients g inside mesh are considered. 


(Va + Va) 2 + (Va + Ve) h2 


V5 2 (h2 + h2) 


(17.11) 


17.3.2 xyLV_RELAX5 P5 B 
hig F Ny 
Relaxation formula is the same as xyLV_RELAX5 P5 A 


(V2 + Vg) h2 + (Va + Vo) he 


V5 = 2 (h2 + h2) 


(17.12) 


17.3.3 xyLV_RELAX5_P5_C 
he =hy=h 


No gradients g inside mesh are considered. 
The formula simplifies, so no g and h terms are necessary. 


_ ve+Vy+Ve+ Ve 
7 4 


Vs (17.13) 


17.3.4 xyLV_RELAX5 P5_D 


The formula also simplifies. 


Relaxation formula is the same as xyLV_RELAX5 P5 _C 


— Vo+V4+ Ve + Ve 


e 4 


(17.14) 
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18 Relaxation formula for node P6 


18.1 Node description 
Right edge of mesh XY. 


18.2 Calculation of relaxation formula 


Laplace equation at node Ps 


Vv? (Vie) py = 


Vow) she Vey) = 
Ox? Oy? 
Pe Pe 


Approximation of partial derivatives of V(,,) at node Ps 


Vox —V Ve—Vs 
(52) Bs a ong a _ geet , Vs— Vo 
Po 


Vo —Ve Vo —V:z 
View me MG e+ Vy-2¥5 
Oy? e iy h2 


Let us substitute approximations to Laplace equation. 


Géo+ , Vs—Ve , Vat Vo 2M _ 
Ng * RBS h2 


Let us find Vg 


Va — Ve V3 + Vo — 2V6 _ G6an+ 


h2 he hex 
Let us multiply both sides 


x 


2,2 
| -Azhe 
We obtain 
V5h2 — Veh + Vghi, + Voh2, — 2Veh2 = —gerthehy 


Vo (2h3, + hy) = (V3 + Vo) he + Vahey + Gor+hohs 
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(18.3) 


(18.4) 


(18.5) 


(18.6) 


(18.7) 


(18.8) 


(18.9) 


(18.10) 


18.3 Final forms of relaxation formula 


18.3.1 xyLV_RELAX5 P6 A 


his hy 
J6a+ # 0 
eee (V3 + Vo) h2 + Vash? + geatheh? (18.11) 
6 a . 
2h2 + h2 
18.3.2 xyLV_RELAX5 P6 B 
hey 
962+ = 0 
(V3 + Vg) h? + Vsh2 
Y= Z (18.12) 
2h2 + h2 
18.3.3 xyLV_RELAX5 P6_C 
(eas eae 
I6x+ # 0 
18.3.4 xyLV_RELAX5 P6 D 
iat 
962+ = 0 
y= Atte (18.14) 
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323 


324 


325 


326 


327 


328 


329 


330 


19 Relaxation formula for node P7 


19.1 Node description 


Left, upper corner of mesh XY. 


19.2 Calculation of relaxation formula 


Laplace equation at node P, 


Vv? (View) p, al 


OV» Vi» 
ou she ead = 
Ox Oy 
Pr P; 


Approximation of partial derivatives of V(,.,) at node P; 


d 


Vs —V; V; —V; Oe 
FPVe,y) ee hg Ve 
7 


2 


= ao 


Viyt+—Vi V7-V: 
(S52) ~ a : Ty _ Va-Vy _ Ity+ 
Oy? o hy he 


Let us substitute approximations to Laplace equation. 


Vga—Vr  g7a— . Va-—Vr . Gry+ 


| + =0 
h2 hy he hy 
Let us find V7 
Ve =? 
Va—Vr | Va—Vr _ Q7a— — Gtyt 
h2 he he hy 
Let us multiply both sides 
272 
| hyhy 


We obtain 


Veh, Vrhi, Vahi, Vehi, = Gta— heh, = G7y+h2hy 


V7 (nz, + hz) = Vih2 + Vgh2 = GTa—hehs 7 Ory phZ hy 
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(19.1) 


(19.2) 


(19.3) 


(19.4) 


(19.5) 


(19.6) 


(19.7) 


(19.8) 


(19.9) 


(19.10) 


332 19.3 Final forms of relaxation formula 
sa: 19.3.1. xyLV_RELAX5 P7_A 
hy £ hy 
Ge—1Gty+ FO 
= Vah?2, + Vahi — g7e—hah2 + gryph2hy 
- (nz + hy) 


(19.11) 


«1 19.3.2 xyLV_RELAX5 P7_B 
he & hy 
G7x—5 J7y+ = 0 


Vihi, + Vahy 
Y= ae (19.12) 


ss 19.3.3 xyLV_RELAX5 P7_C 


hy =hy=h 

Ie-1Ity+ FO 

Va + Ve — gra h 
yew ee (19.13) 
sss 19.3.4 xyLV_RELAX5_P7_D 
hy = hy Sh 
G7x—5 97y+ = 9 
Va + Vi 

V7 = a . (19.14) 
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337 


338 


339 


340 


341 


342 


343 


344 


345 


346 


20 Relaxation formula for node P8 
20.1 Node description 


Upper edge of mesh XY. 


20.2 Calculation of relaxation formula 


Laplace equation at node Ps 


2 
V? (View) », = 9 (20.1) 
PV y) PV y) 
( Ox? Oy? 20:2) 
Ps Ps 
Approximation of partial derivatives of V(,,) at node Ps 
0? Vey) as ee x are _ Va Vo—.2V9 (20.3) 
Ox? 7 he = h2 
Ps 
Vey+—Vs Vg—Vs5 
(Se) _ as 8 iy = Vs — Vg _ IB8y+ (20.4) 
Oy? . hey h2 hy 
Let us substitute approximations to Laplace equation. 
Vi + Vo —2Vg  V5—Ve | Qay+ 
=0 (20.5) 
h2 h2 hy 
Let us find Vg 
Vg =? (20.6) 
V7+Vo—2Vg V5 —Ve J8y+ 
= (20.7) 
h2 he hig 
Let us multiply both sides 
| hehe (20.8) 
We obtain 
Vrhi, + Vohi, — 2Vghi, + Vohz — Veh = —gayrhahy (20.9) 
Ve (hz + 2h2) = (Vz + Vo) he + Vah2 + gayphohy (20.10) 
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347 


348 


349 


350 


20.3 Final forms of relaxation formula 
20.3.1 xyLV_RELAX5 P8 A 


hy & Py 
Jsy+ #9 
ae Vash? + (Vz + Vo) ho + gaythihy 


8 
h2 + 2h? 


20.3.2 xyLV_RELAX5 P8 B 
hy # hy 
J8y+ = 0 


__ V5h2 + (Vz + Vo) hy 
eR 2h? 


20.3.3 xyLV_RELAX5_P8_C 
hy =hy=h 
Qsy+ FO 


_Ve+Ve t+ Vat goysh 
3 


Vg 


20.3.4 xyLV_RELAX5 P8 D 
he =hy=h 
J8y+ = 0 


_v%tVtV 


oe 3 
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(20.11) 


(20.12) 


(20.13) 


(20.14) 
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353 


354 


355 


356 


357 


358 


359 


360 


21 Relaxation formula for node P9 


21.1 Node description 
Right, upper corner of mesh XY. 


21.2 Calculation of relaxation formula 


Laplace equation at node Po 


Vv? (Vey) py =0 


OV» Vi» 
( dL m ead = 
Ox Oy 
Po P9 


Approximation of partial derivatives of V(,,) at node Po 


Ox? 


Or Ves 7 Veen Vs Vo—Va - Vg —Va | 
‘é hy fe 
‘) 


2 


Voy+—Vo Vo—Ve 
FVia,y) ae "hy "Toy : _ Ve ~~ Vo ; G9yt+ 
2 2 
Oy Ps hy hy 


Let us substitute approximations to Laplace equation. 


Vg—Vo  gor+ . Va—Vo . Goyt 


2 a + 


=0 
h2 hy 


Let us find Vo 
Vg =? 
Vg — Vo f Vs—Vo You J9y+ 


hz 2 he hy 
Let us multiply both sides 


x 


242 
| -h hy 
We obtain 


Vehy i Vohi, cs Vehi, o Vohi, = ~99n+hehi, a Joythzhy 


Vo (hz +h?) = Veh? + Veh? + goushahe + goythihy 
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(21.1) 


(21.2) 


(21.3) 


(21.4) 


(21.5) 


(21.6) 


(21.7) 


(21.8) 


(21.9) 


(21.10) 


362 


363 


364 


365 


366 


21.3 Final forms of relaxation formula 
21.3.1 xyLV_RELAX5 P9 A 
hg £ hy 
99x-+599y+ FO 
pe Veh? + Vahi + goothah? + goyrh2hy 


h2 + he 


21.3.2 xyLV_RELAX5 P9 B 
he # Ry 
J9x+99y+ = 9 


Veh2 + Vgh2 
ae 6 Blby 
h2 + h2 


21.3.3 xyLV_RELAX5 P9 C 
ee enn 
J9a-+1I9y+ FO 
_ Ve + Va t+ garth + goyrh 


Yo 2 
21.3.4 xyLV_RELAX5 P9 D 
hy = hy =h 
J9x+ 5 99y+ = 9 
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(21.11) 


(21.12) 


(21.13) 


(21.14) 
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